Mapping the allowed parameter space for decaying dark matter models 
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I consider constraints on a phenomenological decaying dark matter model, in which two weakly 
interacting massive particle species have a small mass splitting, and in which the heavier particle 
decays to the lighter particle and a massless particle on cosmological time scales. The decay param- 
eter space is parameterized by Vk, the speed of the lighter particle in the center-of-mass frame of the 
heavier particle prior to decay, and the decay time r. Since I consider the case in which dark matter 
halos have formed before there has been significant decay, I focus on the effects of decay on already 
formed halos. I show that the «k — t parameter space may be constrained by observed properties 
of dark matter halos. I highlight which set of observations is likely to yield the cleanest constraints 
on i>k — t parameter space, and calculate the constraints in those cases in which the effect of decay 
on the observables can be calculated without iV-body simulations of decaying dark matter. I show 
that for Wk > 5 x 10 3 km s _1 , the 2 = galaxy cluster mass function and halo mass-concentration 
relation constrain r > 40 Gyr, and that precise constraints on r for smaller Wk will require iV-body 
simulations. 

PACS numbers: 



I. INTRODUCTION 

Most of the mass energy of the Universe consists of 
at least two unknown but observationally substantiated 
materials (e.g., Ref. [lj). While there is significant un- 
certainty as to what the "dark energy" that makes up 
3/4 of the Universe might be (for a review, see Refs. 
I ) , there are a number of attractive candidates for the 
~ 1/5 of the Universe made of dark matter. A popular 
class of dark matter candidate is the weakly interact- 
ing massive particle (WIMP), which includes the super- 
symmetric neutralino and the universal extra-dimensions 
Kaluza-Klein photon 0, Hj]. WIMPs are appealing be- 
cause they appear naturally in extensions to the stan- 
dard model (SM) of physics, are thermally produced in 
the early Universe in the quantity required by observa- 
tions, and allow the large-scale structure of the Universe 
to evolve in a manner that appears consistent with ob- 
servations (e.g., Refs. [HIH-El)- 

There are specific predictions on small scales, too. The 
relatively large mass of WIMPs (> 100 GeV) imply that 
thermal decoupling occurs when the WIMPs are nonrela- 
tivistic ( "cold" ) and with a high phase-space density. The 
low speeds in combination with the typically weak scat- 
tering cross sections mean that only the smallest-scale 
perturbations are washed out by the time WIMPs kinet- 
ically decouple from light SM particles [Trj| - [T3j . Hence, 
dark matter halos should be extremely dense at their cen- 
ters and should exist down to ~ Earth mass scales. More- 
over, dark matter density perturbations should evolve as 
a collisionless, pressureless fluid. 

However, it is quite possible that dark matter con- 
sists of something other than standard WIMPs. First, 
while observations are consistent with cold dark matter 



(CDM) on large scales (length scales of > 10's of Mpc), 
these observations do not require a purely cold dark mat- 
ter model. Second, the interpretation of data on smaller 
scales is less clear or is nonexistent (e.g., Refs. flil. [l5T]). 
Third, WIMPs are by no means the only attractive class 
of particle candidates for dark matter. There are a num- 
ber of theories in which the large-scale successes of cold 
dark matter are produced but deviate at the small scales 
where observations either do not exist or are more diffi- 
cult to interpret [16|-[r|. 

In this work, I consider observational constraints on a 
phenomenological decaying dark matter model in which 
a parent particle X experiences two-body decay of the 
form X — » Y + £, where Y is the stable daughter dark 
matter particle and £ is a massless (or nearly massless) 
particle with extremely weak to nonexistent couplings 
with SM particles. The Y particle is only slightly less 
massive than the parent particle X , such that 



M x - M Y ^ , 
M x 



(1) 



For small e, the recoiling particle Y receives a nonrela- 
tivistic speed 



v k /c - 



(2) 
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in the center-of-mass frame of the parent X particle. I 
consider large decay times r > 10 12 s, such that decays 
occur after the first dark matter halos form in the early 
Universe, making this model similar to that considered in 
Sanchez-Salcedo (2003) [20f. Such a model may arise in 
"hidden sector" theories, in which the SM does not inter- 
act with hidden parallel particle sectors (see Refs. [2p - 
[1H for examples of hidden sector theories). Moreover, 
I assume that the parent X particles are massive and 
thermally produced, such that the linear matter power 
spectrum in the early Universe should resemble that of 
CDM. 
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Since this model does not produce SM particles and the 
decay times are long, previously published constraints on 
decaying dark matter do not apply [25l434| . Since the de- 
cay times I consider are long enough that dark matter 
halos should have formed, and formed with CDM-likc 
properties, the main effect of decays is either to inject 
kinetic energy into or eject particles from existing dark 
matter halos. The result of these effects is a change in 
the dark matter density profile in halos and halo mass 
loss. Thus, I consider observational constraints on the 
model from the structure and abundances of dark mat- 
ter halos. In Sec. HH I relate the two free parameters of 
the dark matter model, and r, to the typical speeds 
and time scales of particles in halos, and show how I 
expect halos to respond to decays as a function of the 
decay parameters. In Sec. IIII| I describe the types of ob- 
servations that may constrain this decaying dark matter 
model, and how to translate observations into constraints 
on the decay parameters via the mapping in Sec. [TTJ I 
calculate constraints in those instances in which they can 
be calculated semianalytically. I summarize and discuss 
the results in Sec. HVl 



II. RELATING DECAY PARAMETERS TO 
HALO PARAMETERS 

The phenomenological decaying dark matter model is 
parameterized by just two numbers: the decay time r 
and the kick speed (alternatively, the fractional mass dif- 
ference between dark matter particles) «k- The effects 
of decay on dark matter halos as a function of these pa- 
rameters depends on halo properties, namely the typical 
dynamical time (idyn) and the typical speed of dark mat- 
ter particles in the halo. 

The typical speed of dark matter particles in the halo 
is of order the virial speed, 



virial radius because 



^vir v 



vir 

Rvir/Rvir = COllSt, 



(5) 

(6) 
(7) 



which is the same for all halos. I must select a halo model 
in order to calculate the dynamical time as defined as 
the crossing time at the half-mass radius. Dissipation- 
less cosmological simulations of cold dark matter alone 
(without baryons) show that dark matter halo profiles 
are well described by the Navarro-Frenk- White (NFW) 
profile 



p(r) = 



(S) 



on observable scales (radii from galactic centers r > 
0.01i?vir) [H, [H, H3] • The scale radius r s can be related 
to the i? v ir by 



Rvir / C> 



(9) 



where c is the halo concentration. The concentration is 
expected to be a function of the formation time of the 
halo, and is thus in general a function of the mass of the 
halo [38M40I ] . The scale density p s can be related to virial 
quantities, 



A v p u 



3 ln(l + c) - c/(l + c) ' 



(10) 



If c 3> 1, which is expected for virialized halos, then 
the half-mass radius 



1.65Vc. 



(11) 



The typical dark matter particle speed at such a radius 
is typically 



GM V 



Ry 



(3) 



where the halo virial mass M v i r is defined as the mass 
enclosed within the virial radius i? v j r in which the average 
mass density is equal to the spherical collapse density 
A v p u , 



M v 



4?r 



-A vPu R 



3 

vir ' 



(4) 



Here, p u is the mean mass density in the universe. In a 
ACDM universe, A v (z) = (18tt 2 - 39x - 82x 2 )/tt m (z), 
where x = Q m (z) — 1 and Q m is fraction of the critical 
density p c in the form of matter (p u = Cl m p c ) [351 ] . 

I take the typical dynamical time of a dark matter 
particle in a halo to be the crossing time at the half-mass 
radius of the halo. I do not use the crossing time at the 



Vh 



GM V 



2r h 

(GA v p u ) 1/2 c^ 4 R vil , 



(12) 
(13) 



which is similar to the virial speed. Thus, the typical 
crossing speed at the half-mass radius, and hence, the 
typical dynamical time, is 



tdyn ~ Hi/fh 

« (GA v p u ) 



-1/2-3/4 



(14) 
(15) 



In CDM, a typical galaxy (Af vir ~ 10 12 M Q ) has c « 15, 
yielding v v - 1T = 130 km s _1 and idyn = 500 Myr; a large 
galaxy cluster has M v - lr ~ 1O 15 M0 and c ~ 5, corre- 
sponding to iVir = 1300 km s _1 and t^ ya ~ Gyr; and 
the ultra-faint dwarf galaxies in the Local Group likely 
have virial masses ~ 10 9 M Q (u V i r = 13 km s _1 ) [38l |40T- 
[43] ] . The concentration of those galaxies is unknown, but 
taking c = 30 leads to idyn = 200 Myr. 
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Another interesting time scale is the crossing time at 
the scale radius, t s , since this corresponds to the region 
of the halo at which the profile transitions from p oc 
to p oc r~ 3 . The "scale time" t s is ~ 50 Myr for an 
ultra-faint dwarf galaxy halo, ~ 100 Myr for a typical 
galaxy halo, and ~ 400 Myr for the cluster-sized halo. 
This highlights the point that although I choose one time 
scale to parameterize the dynamical time in a halo, there 
is actually a diversity of dynamical times within a halo. 

Now that I have parameterized the time scales and 
speeds of dark matter particles in halos, I will classify the 
effects of dark matter decays are on halos as a function 
of the decay parameters with respect to the halo param- 
eters. At the present, I will ignore cosmological accretion 
onto halos, and consider the halos to be isolated and in 
equilibrium. 

Case 1 (t > idyn, > v v - lT ) — This is a regime in 
which the decay time is long relative to the dynamical 
time of the halo, and the kick speed is high. The limit 
of relativistic kicks was studied in Refs. |4J[ and [45| in 
other contexts, but even for nonrelativistic kicks greater 
than iv; r , the Y particles are ejected from the system. 
The decays do not directly inject kinetic energy into the 
bound halo, as they will in the following cases, because 
the decay products Y will be unbound to the halo. 

In cases in which the time scale for change in the gravi- 
tational potential is significantly longer than the dynam- 
ical time, particles on regular orbits should conserve adi- 
abatic invariants. This is useful because, if the halo is 
quasi-static, the distribution function (DF) of dark mat- 
ter particles is a function of adiabatic invariants. The 
mass density in the halo thus be calculated using this 
DF if the gravitational potential of the halo is known (cf. 
Ref . [46| ) . In the simplified case of a spherically symmet- 
ric potential with all particles initially on circular orbits, 
the gravitational potential, and hence, the dark matter 
density of the halo may easily be found as a function 
of the fraction of particles that have been ejected from 
the halo. The approximation of dark matter particles on 
circular orbits is not as unrealistic as one might expect; 
cosmological simulations of dark matter halos show that 
there is not a significant phase-space density of highly ra- 
dial orbits, except in the outskirts of the halo [43]. This 
approximation was used in Refs. 44 1 and [45[ , and I show 
analytically how the density profile of dark matter halos 
changes as a function of the fraction of the X particles 
that have decayed, /. 

If the angular momentum of particles is conserved, 
then the initial and final mass distributions are related 



as 



Mi(n)n = M f (rt)rt, 



(16) 



since the specific angular momentum of particles on cir- 
cular orbits is J = (GM (r)r) 1 / 2 , where M{r) is the mass 
enclosed within radius r. Here, i denotes the initial halo 
properties, and f denotes halo properties after a fraction 
/ of X particles have decayed. If the X particles are 
initially on circular orbits, then their orbits do not cross, 



allowing for the following relation: 

M £ (r f ) = (1 - /)Mi(n). 



(17) 



Inserting this equation into Eq. (Ti"6|) . I find that the 
initial and final particle radii are related as 



such that 



= (1 - f)n, 



Mt(r t ) = (1 - /)Mi((l - f)r f ). 



(18) 



(19) 



This changes the mass density of the X particles in the 
following way. The density is 



Pf(rt) 



1 dM f 
Airrf drf 



drf 



(1 - ff dMj 
47rrj dr; 

In the case of an NFW halo [Eq. ((8|)], this implies 

(1 - f) 4 Ps 



Pf(n) = 



(1 ~ f)r t 



(l - f)n 

r s 



(20) 
(21) 
(22) 

(23) 



Thus, the halo retains its NFW form, but with scale ra- 
dius 



r s>f = (l-/) r S)i , 
and with a decreased scale density 
PsJ = (1- /) 4 Ps,i- 



(24) 



(25) 



Thus, if t 3> idyn and ^> v V ir, the shape of the 
dark matter halo will be unchanged, but the scale radius 
will increase, and the mass density and total virial mass 
will decrease. Note that these results are independent of 
«k beyond the fact must be large enough to unbind 
any Y particle from the halo. 

Case 2 (r > tdym < v v - lr ) — In this regime, the halo 
is slowly evolving as a function of time, but the kick 
speeds are small. Unlike case 1 (above), it is difficult to 
calculate analytically the general behavior of the halo. 
Although the decay time scale is long enough that the 
gravitational potential of the halo should evolve fairly 
slowly, the Y daughter particles largely stay within the 
halo, making it difficult to estimate changes to the adia- 
batic invariants or to the gravitational potential. 

However, there are a few general predictions one may 
make. First, since self-gravitating systems have negative 
heat capacity, any slow injection of kinetic energy into 
the halo causes the halo to expand and for the typical 
particle energy to become less negative (see, e.g., Ref. 
46]). Second, there will be some mass loss as some ini- 
tially loosely bound particles will decay to particles that 
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are no longer bound to the halo, and more may be lost 
as the gravitational potential responds to the kinetic en- 
ergy injection resulting from decays. The consequences 
of these effects are to drive down the central halo density, 
the total virial mass, and the typical particle speed. The 
effects will be larger for high and shorter r, and will be 
most pronounced at the halo centers (where the typical 
particle speed is smaller). For small i>k or large r, there 
should hardly be a change to halo properties. 

In general, the behavior of decays in this regime must 
be examined cither by solutions to the Boltzmann and 
Poisson equations or by A-body simulations. 

Case 3 (r < <dyn, «k > v v - 1T ) — This regime, in which 
the decay time is less than the dynamical time, can be 
thought of as being similar to the case of instantaneous 
decay. In this particular case, the kick speed is also quite 
high. In general, most, if not all, of the mass in the 
halo will be ejected. If any mass remains in the halo, 
the system will settle to a new equilibrium within several 
dynamical times, but it is not clear what the structure 
of that halo will be. In the absence of accretion, the 
structure of the halo is fixed for the rest of time. 

Again, quantitative predictions for this regime require 
iV-body simulations. 

Case 4 ( T < idyn, Vk < v v ir) — In this case, the veloc- 
ity perturbations are small and occur in one short epoch. 
This is analogous to the case of high-speed galaxy en- 
counters. And like noncatastrophic high-speed galaxy 
encounters, we can calculate one post-encounter (or post- 
decay) property of the halo: the total kinetic energy in 
the virial radius. If a halo is initially in virial equilibrium, 
the total energy E is related to the kinetic energy K by 
(cf. Ref. 0): 

E = -K. (26) 

If no particles are ejected as a result of the decays, the to- 
tal kinetic energy injected in the halo is AK = M v i t v^/2. 
Thus, the final energy will be 

Et = -K+ AK, (27) 

and after the halo has settled into virial equilibrium again 
after a few dynamical times, the kinetic energy will be 

Kf = -E f (28) 
= K — AK. (29) 

Once again, the negative heat capacity of self-gravitating 
systems means that the net effect of injection of kinetic 
energy into a halo is a decrease in the total kinetic energy 
of the system. 

If Uk is a significant fraction of u v ; r , this approxima- 
tion for the final kinetic energy becomes worse, as decays 
result in mass loss from the system, both directly as in- 
dividual Y particles are created on unbound orbits, and 
indirectly due to the rapidly changing gravitational po- 
tential. 

Even if the final energy of the halo is known, the form 
of the gravitational potential (and hence, the density dis- 
tribution) is not without the help of ./V-body simulations. 




V v vi, 



FIG. 1: Summary of the changes to dark matter halos due to 
late-time decays, as a function of Uk/t>vir and r/td yn . In the 
region marked "CDM-like" , decays have little effect on halos, 
and so halos should resemble those produced by CDM. The 
area marked "catastrophic" is the region of parameter space 
in which halos are almost completely destroyed as a result of 
decays. The middle region labeled "potentially interesting" 
indicates the part of parameter space that is not ruled out 
as catastrophic, nor is it likely to produce halos that exactly 
resemble those predicted by CDM. 



However, just as in case 2, we predict that the halo will 
become less dense, and that the mass loss and density 
changes will be more extreme for larger values of v^. 

There are several points to make about these regimes. 
First, given that galaxies and clusters span a large range 
of halo mass, a fixed r and i>k will place a low mass halo in 
a different regime than a high mass halo. For example, 
if t r~j 1 Gyr and v ~ 100 km s" 1 , ultra-faint dwarf 
galaxies will be in case 1. They will have disintegrated 
by the present since particles will be ejected from the halo 
on time scales far less than the age of the universe, £h ~ 
14 Gyr. However, the halo of an typical (~ 10 12 M Q ) 
galaxy will be in the regime of case 2, and the halo of 
massive cluster will be case 4. Secondly, there will be 
a differentiation in cases throughout a single halo since 
dynamical time scales are much shorter deep within the 
halo. A halo may be in case 2 in the interior, but in case 
4 near the edge of the halo. Last, while one can make 
general predictions for the behavior of halos for decays 
in each regime, quantitative predictions largely require 
iV-body simulations, which we perform and analyze in 
another paper 48]. 

The results are summarized in Fig. [TJ The region 
marked "CDM-like" denotes the r 3> idyn and «k <C v v - 1T 
parts of parameter space in which we expect either al- 
most no decays or almost no kinetic energy injection in 
the halo. For this region of parameter space, the halo 
properties will be nearly indistinguishable from the case 
of no decays. The region marked "catastrophic" has ex- 
tended beyond the parts of parameter space of case 3 
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to lower i>k and higher r, since v\, ~ u v j r can still in- 
duce significant mass loss, and decays on time scales of 
a few dynamical times can affect the halo as it settles 
into its new equilibrium. The region marked "potentially 
interesting" includes bits of cases 1, 2, and 4, and repre- 
sents the part of parameter space whose effects on halo 
structure need to be explored with iV-body simulations 
and constrained with observations. In the next section, 
I show how to relate observational constraints to Fig. 
[U and define regions of t>k — t parameter space allowed 
by observations and which need to be better understood 
with simulations. 



III. OBSERVATIONAL CONSTRAINTS 

There are a number of probes of the distribution of 
matter in the Universe. In this section, I describe a set 
of probes from which constraints on decaying dark matter 
are easiest to infer, and calculate constraints for case 1- 
type decays. We consider observational constraints on 
case 2, 3, and 4 decays in other work [48]. The probes I 
consider are the cluster mass function (Sec. IIII A|) . galaxy 
clustering (Sec. IIII Bp . the existence of small dark matter 
halos (Sec. IIII C|) and the mass-concentration relation 
(Sec. IhTdI . 



A. The Cluster Mass Function 

The cluster mass functions are relatively clean for de- 
caying dark matter studies for both observational and 
theoretical reasons. Individual cluster halo masses are 
substantially easier to determine than the halos of indi- 
vidual galaxies, and the assignment of a cluster of galax- 
ies to a single dark matter halo is unambiguous. More- 
over, since most of the energy in decays goes to nonrel- 
ativistic massive particles, the background evolution of 
the Universe is unchanged. Hence, the growth function, 
which depends only on redshift and the background cos- 
mology in a ACDM universe, is unchanged relative to 
ACDM predictions (see Refs. [H, [50| for the relativis- 
ts case, in which the background evolution does change) 
(5lj . Thus, main change to the cluster mass function rel- 
ative to that predicted from ACDM models occur as a 
result mass loss from the halos due to decay. This makes 
the mass functions simple to interpret with respect to 
ACDM predictions. 

There are two different ways to use the cluster mass 
function to constrain the decaying dark matter param- 
eter space. I will show how the z = mass function 
(see Refs. [j| l52l - [5t| for recent observations and mass- 
function calculations) constrains the decay time r in the 
case 1 virial regime. This method is easy to generalize 
to other parts of decay parameter space. Second, I show 
how upcoming cluster surveys (e.g., Sunyaev-Zel'dovich 
observations in the microwave bands [57]) that are sensi- 
tive to the redshift evolution of the cluster mass function 



can be used to constrain decays. 



1. The local cluster mass function 

Cluster mass functions derived from X-ray and optical 
surveys of the local Universe have been used to constrain 
Q m and as , the rms amplitude of fluctuations in the lin- 
ear density field at z = if smoothed on 8 h~ x Mpc 
scales, where h is the Hubble parameter 0, IH^l ■ The 
constraints on VL m and as from the local cluster mass 
function are consistent with those found with other cos- 
mological probes, including the cosmic microwave back- 
ground (CMB) [58], the Lyman-alpha forest B, w eak- 
lensing power spectra, and galaxy clustering [8j, |59M65| . 

The consistency of the constraints on f2 m and as from 
different epochs constrains decaying dark matter for the 
following reason. The standard lore of halo formation 
is that halos form when the amplitude of the density 
perturbation (alternatively, the amplitude of the linear 
matter power spectrum) smoothed on a distance scale 
R (or alternatively, a mass scale M = 47T/9 u i? 3 /3), ex- 
ceeds the overdensity required for halo collapse [66|, [67} ■ 
The halo number density dn/dM is then related to the 
volume of space in which the overdensity exceeds the col- 
lapse density. If the smoothed amplitude of fluctuations 
is small, then there are few virialized halos of that given 
mass; if the amplitude is large, then there are many ha- 
los of a given mass. The amplitude of the linear matter 
power spectrum is controlled by the initial amplitude of 
fluctuations from inflation as well as f2 m , but is usually 
parameterized in terms of as and Q m . 

The number density of clusters, the rarest of virial- 
ized objects, is quite sensitive to as and f2 m , with more 
clusters expected in a high-erg, high-Sl m cosmology. A 
~ 20% change in as can yield order-unity changes in 
the cluster halo number density (e.g., Ref. [68]). De- 
cays cause mass loss in halos, thus reducing the number 
density of halos above fixed mass. If one were to infer 
ACDM parameters from the local cluster mass function 
if dark matter were decaying, one would infer artificially 
small values of as and f2 m relative to what one would 
infer from probes of earlier epochs (e.g., the CMB and 
Lyman-alpha forest). The fact that as and Q m inferred 
from the local cluster mass function are similar to those 
inferred from other probes limits decay parameter space, 
although the error bars for the cluster mass functions 
(and hence, those for the inferred values of as and f2 m ) 
are quite large [H, [H, [H| . 

One can estimate conservative limits on decay parame- 
ter space using the following method. Since decays lower 
the inferred as and Q, m of clusters relative to those in- 
ferred from the CMB, I find which combination of and 
r would cause the 2-a upper limits on O m and as from 
the CMB to yield a z — cluster mass function that is 
barely consistent to the 2-a lower limits on those param- 
eters. In the instance of case 1 of Sec. HH this can be 
done analytically. Here, I show the limits one can set on 
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FIG. 2: Halo number density n(> M) as a function of mini- 
mum halo mass M at z — 0. The dotted (blue) line represents 
the number density of halos in a cosmology with Q m = 0.318 
and ag — 0.868 (with both quantities being 2-a above the 
WMAP-5 mean values), the solid (green) line represents the 
halo number density for a WMAP-5 cosmology (fl m = 0.258, 
as = 0.796), and the dashed (red) line shows the number den- 
sity of halos in a cosmology with Q m (=0.198) and as (=0.724) 
2-a below the WMAP-5 mean values. These number densi- 
ties were generated using the Tinker et al. [6|| mass-function 
fitting formulas. The dot-dashed lines show halo number den- 
sities (top to bottom) once a fraction / = 0.1, 0.2, and 0.3 of 
the X particles have decayed to fast Y particles, assuming a 
high-f2 m , high-(78 cosmology. 



r assuming that is several times the typical cluster 
virial speed (« v ; r ~ 10 3 km s _1 ). 

I take the Vt m and as and the associated error bars from 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
five-year data set as the baseline, since this is the ear- 
liest epoch from which f2 TO and the amplitude of the 
matter power spectrum can be inferred [58|. The cen- 
tral values and 1-a error bars from the WMAP-5 6- 
parameter ACDM fits are Q m = 0.258 ± 0.030 and 
o-a = 0.796 ± 0.036. I use the Tinker et al. (2008) [H 
mass function to estimate the comoving number density 
n(> M) of halos above a mass threshold M at z = 
for fl m = 0.318 and a 8 = 0.868 (2-a above the WMAP- 
5 mean values) and f2 m = 0.198 and as = 0.724 (2-a 
below), keeping the Hubble parameter h, the primordial 
slope of the matter power spectrum n s , and the baryon 
fraction f2f,/i 2 fixed to the WMAP-5 central values, as- 
suming a flat ACDM cosmology. The cluster mass func- 
tion is far less sensitive to those parameters than to £l m 
and as- The mass functions are shown in Fig. [5J 

Then, assuming there is a one-to-one mapping between 
a ACDM halo virial mass M; (calculated using the spher- 



ical top-hat overdensity described in Sec. |TTJ) for the high- 
f2 m , high-erg cosmology and the virial mass Mf after a 
fraction / of the X particles in the halo have decayed, 
I map the high-f2 m , high-erg ACDM comoving number 
density n\ to the comoving number density of halos after 
a fraction / of X particles have decayed 



Uf(> M { ) 



,dni(Mi(M^) dMi 



dMi 



dM[ ' 



(30) 



The resulting comoving number densities are shown as 
the dashed lines in Fig. [2]for / = 0.1, 0.2 and 0.3. 

The comoving number density of cluster halos with 
M vir > 1O 14 /i" 1 M (roughly the mass threshold for clus- 
ter studies) if / = 0.3, corresponding to r ps 40 Gyr, is 
nearly identical to that of the low-O m , low-erg cosmology, 
although the mass function has a steeper slope at larger 
masses ~ 5 x lO 14 /i _1 Af . For slightly smaller values of 
/ (ps 0.25), the decay mass function lies above the low- 
f2 m , low-erg mass function for M v ; r < 10 15 h~ 1 MQ. Thus, 
the upper limit on / is near 0.3, such that r > 40 Gyr 
for i>k > a few times the cluster virial speed of the most 
massive observed clusters. Since clusters are observed 
to have M v i r < 2 x 10 15 Af Q , this calculation applies to 

> 5000 km s~\ 

If the decay parameters relative to the cluster dynam- 
ical parameters are in any case other than case 1 of Sec. 
III1 one will need to perform iV-body simulations of decay 
in halos in order to map decay parameters to the halo 
mass loss. 

The / — 0.3, t ~ 40 Gyr lower limit for case 1 is likely 
to be too conservative, since I have assumed a one-to- 
one mapping between a ACDM halo of mass Mj and a 
halo in the decay case with mass Mf. However, to do 
this sort of analysis self-consistently, one should perform 
cosmological TV-body simulations to take into account the 
streaming of dark matter particles. 



2. Redshift evolution of the mass function 

The evolution of the cluster mass function is a promis- 
ing way in which to explore evolution in the dark-energy 
equation of state w, since the evolution in the comoving 
number density of cluster mass halos is quite sensitive 
to the growth function and to the Hubble constant as 
a function of redshift, H(z). While current constraints 
on w from cluster mass functions are weak Q, ongoing 
and future Sunyaev-Zel'dovich surveys and follow-up ob- 
servations should provide much better constraints (e.g., 

Ref. HH!). 

The same reasons why the evolution of the cluster mass 
function is a useful probe of dark energy also make it 
a useful probe of deca ys, which was originally pointed 
out in Refs. [45[ and [43|. In the right-hand panel of 
Fig. [31 I show the z = comoving number densities of 
halos for the high-fi m , high-ag; mean WMAP-5; and low- 
Q m , low-erg cosmologies, as well as the comoving number 
density of halos if t = 50 Gyr in case 1 of Sec. HH fl m = 
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FIG. 3: Evolution of the number density of massive halos as 
a function of time. The lines have the same meaning as in 
Fig. [2] and the dot-dashed line represents a high-f2 m , high- 
erg cosmological model with > « v ir of clusters and r = 50 
Gyr. 



0.318, and erg = 0.868. The middle and left panels show 
the comoving number densities at z — 0.5 and z = 1, 
respectively. One can see that the z = mass function 
for r — 50 Gyr lies between the low-f2 m , low-erg and mean 
WMAP-5 mass functions, but it lies between the mean 
WMAP-5 and high-f2 m , high-erg cosmologies by z = 1. 
The comoving number density of halos above lQ 1A h~ l MQ 
in the r = 50 Gyr cosmology only increases by a factor 
of 5 since z = 1, even though it grows by a factor of 
~ 8 in the high-fi m , high-erg cosmology, a factor of ~ 11 
in the mean WMAP-5 cosmology, and a factor of 16 in 
the low-f2 m , low-erg cosmology in the absence of decays. 
Thus, a signature of decaying dark matter would be a 
slower growth of structure with redshift than expected 
based on the z = mass function. 

A caveat is that a non-cosmological constant (A) model 
for dark energy could also change the evolution of n(> 
M) above a fixed mass threshold. However, it may be 
possible to break the degeneracy between the two by 
looking at the evolution of the shape of the cluster mass 
function. 



B. Galaxy clustering 

The clustering of galaxies is often used to constrain 
cosmological parameters. On large, linear scales (which, 
at z = 0, corresponds to comoving wavenumbers k < 
O.lh Mpc -1 ), clustering is usually expressed in terms of 
the galaxy power spectrum P gg (fc), the Fourier transform 
of the galaxy autocorrelation function. On smaller, non- 
linear scales, the clustering is usually expressed in terms 
of the real-space galaxy two-point correlation function 
£(r). Galaxies are biased tracers of the matter distribu- 
tion, but the matter power spectrum P(k) (or correlation 
function) can be extracted from galaxy clustering using 
analytic and empirical models of how galaxies populate 
dark matter halos. 

In this section, I describe how P(k) and £(r) could be 



used to constrain decaying dark matter models, but a 
detailed calculation is far beyond the scope of this pa- 
per. I will give only a crude, order-of-magnitude limit 
on decaying dark matter parameter space based on the 
measured £(r) from the Sloan Digital Sky Survey (SDSS) 
[70| . A proper analysis of the limits on decaying dark 
matter would require solving the Boltzmann and Einstein 
equations for the linear matter power spectrum and the 
collapse criteria for halo formation, as well as A^-body 
simulations to better explore the nonlinear regime. 

The linear matter power spectrum can be determined 
from first principles by solving the Einstein equations 
for first-order perturbations from the background met- 
ric and the Boltzmann equations governing the evolution 
of the distribution functions of all the contents of the 
Universe. Calculations typically include CDM, baryons, 
massive neutrinos, and a cosmological constant A for the 
dark energy. Relative to a ACDM + baryon power spec- 
trum, the ACDM + baryon + neutrino power spectrum 
is suppressed on scales corresponding the the neutrino 
free-streaming scale [71| . This is because neutrinos can 
escape the matter density perturbations if their typical 
speed is larger than the escape speed from a perturba- 
tion, which reduces the density perturbation. In decaying 
dark matter models, the X particles continuously source 
a free-streaming ^-particle populations. Thus, the den- 
sity perturbations on small scales (below the typical free- 
streaming scales of the Y particles) should be damped 
relative to the ACDM expectation. However, the power 
spectrum on large scales should be unaffected by this 
damping. 

The scale-dependence of the damping of the power 
spectrum due to decays is what may make the effects of 
decays distinguishable from non-A dark-energy models. 
Unless dark energy clusters, the primary effect of a non- 
A dark energy is to alter the growth function, which is 
scale independent. Thus, to distinguish between varying- 
w and decaying dark matter models, one will need to de- 
termine the evolution of the shape of the linear matter 
power spectrum (or the cluster mass function). 

On nonlinear scales, the two-point galaxy correlation 
function £(r) is often used to quantify galaxy clustering, 
instead of the Fourier-space analog preferred on larger 
scales. Since galaxies are biased tracers of the underly- 
ing matter field, considerable effort has gone into under- 
standing the connection between different galaxy pop- 
ulations and the matter field. This is also relevant on 
larger scales, to map between the galaxy power spec- 
trum and the linear matter power spectrum. The cor- 
relation function is often analyzed in the context of 
halo-occupation distribution (HOD) models, analytic but 
simulation-calibrated relations describing the distribu- 
tion P(N\M) of the number N(M) of a galaxies of a 
certain type contained within a halo of mass M (z^-tzl] ■ 
P(N\M) is determined empirically, along with cosmolog- 
ical parameters, for the set of galaxies being analyzed. 
There are some degeneracies between the HOD and cos- 
mological parameters, but these can be broken by com- 
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bining the £(r) analysis with different galaxy statistics of 
the same set of galaxies, or different observations alto- 
gether [zllzl. 

HODs make the interpretation of £(r) simpler, so I can 
show why it is difficult to pull easy constraints on decay- 
ing dark matter models from the correlation function, 
and to show what would be necessary in order to use 
£(r) to quantitatively constrain the dark matter model. 
The following description closely follows that of Ref . [72J ■ 

The galaxy correlation function can be broken into 
parts, 



£M = l + £ih(r)+6h(r), 



(31) 



where ^lh^) is one- halo term, the correlation function of 
galaxies within the same halo, and £,2h(r) two- halo term, 
the correlation function of galaxies in different halos. The 
one-halo term is 



l + 6h(r) 



1 



27rr 2 n 2 



dM 



aln (N(N-l)) M 
dM 2 



F' 



2R vix (M) \2R 



(32) 



where n g is the number density of galaxies, dn/dM is the 
halo mass function, (N(N— 1)) m/2 is the average number 
of pairs of galaxies of a specific type within a halo of mass 
M, and F' (r /2R v - lr ) is radial distribution of the galaxy 
pairs. n g is a measured quantity and (N(N — 1))m/2 
is inferred in the parameter estimation of the correlation 
function. However, dn/dM and F 1 (r /2R v - a ) are taken 
from dissipationless iV-body simulations. 

The two-halo term, to lowest order, is given by 



6h(r) 



1 

2^2 



30 dkP ^ {k) e h 2^L. 

56 fc r 



(33) 



where 



P^ = P m (k)- 



dM—(N(M))b h (M)y g (k,M) 



(34) 



Here, P gg is the two-halo galaxy power spectrum, and 
P nl is the nonlinear matter power spectrum, which is 
taken from iV-body simulations. (N(M)), the average 
number of galaxies in a halo, is determined empirically, 
but 6h, the bias of the dark matter halo power spectrum 
with respect to the matter power spectrum, and y g {M, k), 
the Fourier transform of the galaxy density profile in ha- 
los, are taken from iV-body simulations. 

The problem with trying to find precision constraints 
on the decay parameters from £(r) is that the interpre- 
tation of thereof requires detailed iV-body simulations 
to determine P n \k), dn/dM, F'(r/2i? vir ), b g (M), and 
y g (M, k), although a few of those quantities may be found 
empirically with weak lensing (e.g., Refs [7^,|8Q|). These 



simulations are generally only performed for "standard" 
cosmologies, and it is not clear how dark matter decays 
will affect the power spectrum or correlation functions, 
or interpretation thereof. Simulations using ACDM cos- 
mologies show that cosmological parameters are some- 
what degenerate with the halo occupation models for the 
galaxies, but the degeneracies have only been explored in 
a limited set of cosmological models 0, [7(| . 

However, there is general consistency of the fits to the 
correlation function of local SDSS galaxies with ACDM 
models down to halo masses ~ 10 12 M Q [70j . This sug- 
gests that decaying dark matter parameters could be 
well-constrained using the galaxy correlation function. 
In the absence of A^-body simulations to calibrate HOD 
models as a function of decaying dark matter parame- 
ter space, I constrain the parameter space by requiring 
that ~ 1O 12 M ACDM halos lose less than half their 
mass due to decays. This roughly restricts r > 30 Gyr 
for > v v i r pa 130 km s^ 1 , which are case-l-type con- 
straints. It should be noted that this is a highly approxi- 
mate constraint; more precise constraints require cosmo- 
logical A^-body models to calibrate HODs in the case of 
decaying dark matter. 

In addition to simulations to constrain HOD models for 
decaying dark matter cosmologies, one may constrain the 
decay parameter space observationally with gravitational 
lensing. From gravitational lensing, one may determine 
the typical halo mass for a specific type of galaxy. This is 
an additional constraint to the HOD. Main halo masses 
have been found for galaxies binned by stellar mass in 
the SDSS [791 ] ■ and should be even better-determined in 
upcoming large all-sky surveys (e.g., the LSST (8l|). 



C. Existence of Small (~ 1O 9 M ) Halos 

The smallest virialized halos for which there is obser- 
vational evidence of existence arc of order 10 9 M Q . Halos 
of such small size are observed in two different ways, and 
in both cases, the halos are actually subhalos. Strigari 
et al. fi| find that all faint dwarf galaxies in the Milky 
Way halo have mass interior to 300 pc of ~ 10 7 M Q re- 
gardless of luminosity. This corresponds to a virial mass 
of ~ 10 9 Mq if the mass profile is extrapolated beyond 
the stellar component. 

Vegetti et al. 82] find evidence for a subhalo of mass 
M sub = (3.51 ±0.15) x 1O 9 M in the double Einstein ring 
system SDSSJ0946+1106. This is the mass within the 
tidal radius, so the subhalo was likely somewhat larger 
before it fell into the larger halo of the elliptical galaxy, 
the lens system. 

The existence of these small halos can set limits on a 
broad swath of decay parameter space. The most strin- 
gent limit is set if one requires that the number density 
of 10 9 MQ-mass halos be the same as predicted by CDM, 
and that the density profiles within those halos not be sig- 
nificantly disturbed due to decays. Recalling from Sec. 
[IT] that the virial speed of such halos is ~ 13 km s _1 , one 
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can roughly exclude the decay parameter space above 
«k ~ 10 km s^ 1 and r < a few times in- However, there 
are no z = observational constraints on the low-mass 
halo mass function. 

As there is no measurement of the number density or 
correlation function of halos below ~ 10 12 M Q , there is 
considerably more freedom in the decay parameter space. 
For example, one could imagine that the 1O 9 M0 halos 
that have been observed at z = are the remnants of 
larger halos that have been heavily disturbed by decays. 
If today's ~ 10 9 M© halos are parented by halos that were 
initially 10 M , then parent halos would have initially 
had u v ir ~ 30 km s _1 and idyn ~ a few hundred Myr. 
For such halos to lose ~ 90% of their mass, > w V ir and 
t < tn. This opens up much more allowed decay pa- 
rameter space than the more previous, more conservative 
limit. One could imagine halos of up to ~ 10 11 M Q de- 
caying to ~ 10 9 M Q halos; in this case, short decay times 
(t < tii) and large kick speeds (vk > u V ir ~ 60 km s _1 ) 
are required. Pushing against this limit is the fact that 
~ 10 M© halos appear to have correlation functions 
consistent with ACDM cosmology. The virial speed of 
1O 11 M halos is a significant fraction of the virial speed 
of 10 12 Mq halos, so such halos would likely be quite dis- 
turbed, the degree of which can only be ascertained with 
simulations of the cases 2 and 4 of Sec. [IT] 

In summary, the existence of ~ 10 9 M Q halos excludes 
the parameter space «k > 60 km s _1 and r < in- Stricter 
constraints come from requiring the mass function of 
~ 10 9 — lO n M0 halos to resemble that which is predicted 
from ACDM models, although I emphasize that there are 
no observations that require these stricter constraints. 
These constraints are similar to those from galaxy cor- 
relation functions (Sec. IIIIB[) . For < 1 km s _1 or 
t > 10£h, halos with masses > 10 9 A/ Q should look and 
cluster like CDM halos. 



D. Mass-Concentration Relation 

ACDM predicts a relationship between the concentra- 
tion c (Eq. ©) and halo mass, in addition to the shape of 
the dark matter density profile in halos. Low-mass halos 
are expected to be more concentrated than high-mass ha- 
los, since smaller halos form earlier when the Universe is 
more dense. Concentrations are expected to be higher at 
the present than at higher redshift for fixed halo mass. 
These trends have been found in a number of A-body 
simulations (e.g., Refs. J1J1 |H |H, HJH) and ob- 
servations (e.g., Refs. [86- §3), although there is some 
disagreement in the details. 

As with the cluster mass function in Sec. IIII Al I 
can place constraints on decaying dark matter models 
by considering the range of mass-concentration relations 
allowed in WMAP-5 cosmology. High-f2 m , high-erg cos- 
mologies produce higher concentrations for fixed mass 
than low-f2 m , low-erg cosmologies, since halos collapse 
earlier for fixed mass if the amplitude of the matter 
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FIG. 4: Halo mass-concentration relations at z = 0. The 
solid lines are predictions from the JV-body simulations of 
Maccio et al. @ for WMAP-1 (upper) and WMAP-3 (lower) 
cosmologies. The short-dashed and dotted lines are mass- 
concentration relations observed in galaxy clusters [88| and in 
galaxy groups and clusters [83 ■ The symbol in the upper- left 
corner represents the intrinsic scatter in the simulations and 
l-<7 errors in the mean observed mass-concentration relations. 
The shaded region shows the mass-concentration relation in 
Ref. [89( . The long-dashed lines show the mass-concentration 
relation for high-fk dark matter models after (top to bottom) 
a fraction / = 0.1, 0.2 and 0.3 of the X particles have decayed. 



power spectrum is higher. In Fig. 2J I show the mean 
mass-concentration relation for WMAP-1 (upper solid 
line; Q m = 0.299, er 8 = 0.9) and WMAP-3 (lower solid 
line; Q, m = 0.238, erg = 0.75) cosmologies from the dark 
matter-only simulations of Maccio et al. [43[ . The values 
of f2 m and erg bracket the 2-cr values of fl m and erg of 
the mean WMAP-5 cosmology, and the WMAP-5 mass- 
concentration relation should lie between the WMAP-1 
and WMAP-3 relations [H, H| EEO] ■ The error bar in the 
upper- left corner of Fig. [4] shows the intrinsic scatter in 
the relation, which is large. 

As in Sec. IIII Al I set conservative constraints on case 1 
virial parameters (t>k > v v - 1T , r > idyn, where the halo pa- 
rameters correspond to the most massive halos observed) , 
by considering the effects of decay on the observables by 
assuming a high-f2 TO , high-erg (WMAP-1) cosmology. I 
use Eqs. (j2~4"|) and (|2"5"]) to determine the mass and con- 
centration of a halos after a fraction / of A-particles 
have decayed. The new mass-concentration relations are 
shown in Fig. @] with long-dashed lines, with (top-to- 
bottom) / = 0.1,0.2 and 0.3. For / > 0.2 (r > 60 
Gyr), the mean mass-concentration relation lies below 
the WMAP-3 relation from Ref. [Hj|. 

Unlike Sec. IIII Al I use the observed mass- 
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concentration relation to constrain the decay parame- 
ter space. I show mass-concentration relations inferred 
from observations in Fig. [4] The short-dashed line shows 
the mass-concentration relation found by Ref. [88] using 
published masses and concentrations of galaxy clusters. 
The masses and concentrations of these clusters were de- 
termined using a variety of different observations (e.g., 
strong lensing, weak lensing, cluster galaxy dynamics, 
X-ray temperature profiles). The dotted line shows the 
mass-concentration relation inferred from X-ray temper- 
ature profiles of elliptical galaxies, galaxy groups, and 
galaxy clusters (87| . The hatched region shows the \-a 
range of mass-concentration relations from weak lensing 

Of the three observational data sets, the weak lens- 
ing data set should be least affected by selection effects. 
It has been shown in simulations that strong-lens sys- 
tems are biased towards high concentrations for fixed 
mass relative to the population of halos as a whole (92i — 
|94| . Comerford and Natarajan [88[ show that the mass- 
concentration relation of strong-lens systems in the sim- 
ulations of Hennawi et al. [Hf is a close match to their 
mass-concentration relation from clusters. Halo mass 
profiles can only be reasonably determined from X-ray 
data if the halo is relaxed; relaxed halos have higher con- 
centration than halos as a whole [40]. However, selection 
biases for the X-ray probes of halo properties have not 
been quantified. 

Thus, I compare the mass-concentration relation from 
the decaying dark matter models to the weak-lensing re- 
lation. I find that for / = 0.3 (t = 40 Gyr), the decaying 
dark matter mass-concentration relation lies well below 
the mean relation from weak lensing 89]. This sets the 
lower limit on the allowed value of r for high t^, where 
t>k > iVir f° r the largest clusters probed in the observa- 
tions (tvj r > a few xlOOO km s _1 ). This limit is nearly 
identical to that obtained from the cluster mass function 
in Sec. fOAl 



IV. DISCUSSION 
A. Summary 

The constraints are summarized in Fig. [5j At near- 
relativistic to relativistic v^, the only allowed region of 
parameter space is r > 120 Gyr, which is nearly indistin- 
guishable from CDM, as only ~10% of X particles will 
have decayed by the present. The region with r < 20 
Gyr yields catastrophic destruction of halos on all scales, 
while 20 Gyr < t < 120 Gyr is ruled out by the shape 
of the CMB temperature power spectrum j95| , although 
additional constraints may be possible using luminosity 
distances of supernovae [961 ]. 

For kick speeds of order > 5000 km s _1 , both the clus- 
ter mass functions and the mass-concentration relation 
restricts t > 40 Gyr. A calculation of the cosmolog- 
ical Boltzmann and Einstein equations including decay 




v k [km s" 1 ] 

FIG. 5: Diagram of allowed and excluded «k — t parameter 
space. See text for descriptions. 



to find the linear matter power spectrum, and iV-body 
simulations to explore the nonlinear regime, should yield 
even better constraints on r in this regime. 

For 60 kms _1 <u k < 5000 km s~\ r > 30 Gyr based 
on £(r) and the existence of ~ 10 9 M© halos. The ex- 
cluded region is marked (1) in Fig. [5] Better constraints 
on r for this range of i\ using £(r) will come from cos- 
mological simulations of structure formation that include 
decay. It may be easier to constrain this part of de- 
cay parameter space using mass-concentration relations 
or the cluster mass function. In this case, one can do 
simulations of isolated halos in the regimes r < td yn 
and Vk < fvir & t > t^ yn to map the ACDM mass- 
concentration relation and cluster mass function to those 
with cosmologies with decaying dark matter. We perform 
such simulations and will present our findings in another 
paper [48j . 

Any r is allowed for 1 km s _1 < < 60 km s _1 unless 
one uses the more conservative constraints on the param- 
eter space based on the existence of a few ~ 10 9 M Q halos. 
The additional exclusion zones are marked (2) & (3) in 
Fig. [5l corresponding to allowing 10 9 M Q halos to result 
from decays of halos that initially had mass 10 10 M Q and 
from requiring that the halo number density n(> 1O 9 M0) 
be nearly indistinguishable from ACDM, respectively. I 
emphasize that there are no observations that definitively 
exclude these additional regions of parameter space. In 
general, these regions of parameter space can be bet- 
ter constrained theoretically by the aforementioned cos- 
mological simulations, in order to relate the observed 
£(r) to decaying dark matter parameter space. It may 
be possible to constrain some of this parameter space 
using the mass-concentration relation, as weak lensing 
has already probed this relation down to halo masses 
Af v ; r < 10 12 M Q . Future deep surveys covering a large 
fraction of the sky (e.g, LSST [8l|) should improve and 
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extend measurements of the mass-concentration relation 
to lower halo masses. 

For % < 1 km s" 1 , there will be negligible deviations 
from ACDM cosmology on ~ 10 9 M Q halo scales and up, 
and negligible differences on all scales if t > 100 Gyr, 
for the simple reason that structure does not look signifi- 
cantly different from ACDM if almost no particles decay 
or if the daughter particles receive a kick small relative to 
the virial speed. This part of parameter space is labeled 
"CDM-like" in Fig. [5] for those regions. If halos smaller 
than 10 9 M Q are observed, this would push the "CDM- 
like" line at Vk < 1 km s _1 down to smaller v^. The 
discovery of smaller (sub)halos is likely in future obser- 
vations of strongly lensed galaxies, since flux anomalies 
and time delays are sensitive to the su bhalo mass function 
down to near-stellar mass scales 97l - ll0l| . The absence 
of smaller subhalos would have profound implications for 
the nature of dark matter. 



B. Why I am not using the dark matter density 
profile as a constraint 

Dissipationless A^body simulations of structure forma- 
tion suggest that the density profiles of dark matter halos 
should be well described by the NFW profile of Eq. © 
on any observable scale. Although decays in the case- 
1 regime to not alter the structure of dark matter halos, 
decays in other regimes generically change the halo struc- 
ture due to kinetic-energy injection. Thus, in principle, 
the density profile in dark matter halos would appear to 
be a good test of the decaying dark matter model. 

However, observations of galaxies and clusters, from 
the smallest ultrafaint dwarfs in the Local Group to the 
largest virialized halos in the local Universe, show signif- 
icant scatter in the density profile on scales r < r s . Some 
observations indicate cores in the density profile, others 
show cusps shallower than NFW, some show that the 
NFW profile is a reasonable fit, and others show profiles 
that are more sharp ly cusped than NFW (for a summary , 
see Refs. (ToMosj on Milky Way dwarfs, 0, [Toi, fT07j 
on the ro tation curves of low-surface-brightness galaxies, 
1081 - 11 o f on the rotation curves of massive spiral galax- 



ies, 



on lensing an d kinematics of elliptical 
galaxies, and [80l M. \umm | on various probes of the 
dark matter distribution in galaxy groups and clusters). 

While there are likely to be some observational system- 
atics that affect the dark matter profile fits, perhaps the 
most significant issue in interpreting these data is that 
there are serious theoretical systematics. The NFW pro- 
file emerges in simulations of structure formation without 
baryons. However, galaxies dominate the gravitational 
potential in the inner parts of the dark matter halos (with 
the possible exceptions of the ultrafaint dwarfs in the 
Local Group and some low-surface-brightness galaxies). 
There is evidence from simulations that the dark mat- 
ter responds to various processes associated with galaxy 
evolution (feedback, star formation rate, gas cooling, 
smooth or clumpy baryonic infall), but it is not clear 



which processes dominate [I22l 4l28]. Dark matter ha- 
los in A-body simulations with CDM and baryons range 
from ha ving cored density profiles [129l | to quite cusped 
profiles |l30l Il3l| . At this point, the interpretation of 
simulations is descriptive rather than predictive. What 
is needed in order to interpret observations are concrete 
theoretical predictions, not only for CDM, but for other 
types of dark matter. 

Thus, even though dark matter decay is likely to 
change the density profile of dark matter in halos from 
the form of Eq. ©, the observed density profiles will 
not constrain decaying dark matter parameter space un- 
til there are firmer theoretical predictions for the dark 
matter distribution in halos in the presence of a bary- 
onic galaxy. This is especially important because decays 
affect the central regions of halos, the home of galaxies, 
more strongly than the outside regions for fixed r and 
i>k, since the dynamical time scales and typical particle 
velocities are lower. 



C. Future directions 

The easiest way to constrain more of the decay param- 
eter space is with the mass-concentration relation and 
cluster mass functions. Future observational data will 
improve constraints on the decay parameter space. Up - 
coming deep optic al all-sky surveys such as DES |l32j |. 
PanSTARRS [l33| . and LSST [81( will allow for a bet- 
ter determination of the mass-concentration relation via 
weak lensing. The redshift-dependent cluster mass func- 
tion may be determined with next-generation optical 
all-sky sur veys, as wel l as Sunyaev-Zel'dovich surveys 
[SPT PH; ACT [HI) and X-ray surveys (eROSITA 
|l36j |). These data sets, along with CMB measurements 
with Planck, will provide better constraints on cosmo- 
logical parameters, which will further constrain the de- 
cay parameter space. The effects of decay on the mass- 
concentration relation and cluster mass functions is rel- 
atively easy to quantify since these constraints do not 
require cosmological simulations. In this work, I showed 
how to constrain r in the case that is significantly 
larger than the virial speed of the largest galaxy cluster. 
To constrain other parts of the decay parameter space, 
one may do noncosmological simulations of isolated dark 
matter halos. This will allow one to determine the change 
to the halo structure and mass loss as a function of decay 
parameters, and is investigated in Ref. (48[. 
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